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Abstract 

In string or M theories, the spontaneous breaking of 10D or 11D Lorentz symmetry 
is required to describe our space-time. A direct approach to this issue is provided 
by the IIB matrix model. We study its 4D version, which corresponds to the zero 
volume limit of 4D super SU (N) Yang- Mills theory. Based on the moment of inertia 
as a criterion, spontaneous symmetry breaking (SSB) seems to occur, so that only one 
extended direction remains, as first observed by Bialas, Burda et al. However, using 
Wilson loops as probes of space-time we do not observe any sign of SSB in Monte Carlo 
simulations where N is as large as 48. This agrees with an earlier observation that the 
phase of the fermionic integral, which is absent in the 4D model, should play a crucial 
role if SSB of Lorentz symmetry really occurs in the 10D IIB matrix model. 
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1 Introduction 



Matrix models 0, |] are considered the most promising candidate for a nonperturbative 
definition of string/M theories. They may play an analogous role as lattice gauge theory 
does in quantum field theory. One of the most fundamental questions that can be addressed 
using these models is the issue of spontaneous breakdown of Lorentz invariance, which is 
required to occur in order for these theories in 10 (or 11) dimensions to describe our four- 
dimensional space-time. For early works which address this issue using string field theory, 
see Refs. @. 

The IIB matrix model (or IKKT model) 0, which is conjectured to describe type IIB 
superstrings nonperturbatively, is a supersymmetric matrix model composed of 10 bosonic 
matrices and 16 fermionic matrices, which can be obtained formally by taking the zero- 
volume limit of 10D super SU(iV) Yang-Mills theory. [] This model is particularly suitable for 
the study of spontaneous symmetry breaking (SSB) of Lorentz invariance,^ since it is mani- 
festly invariant under SO(10) transformations, which transform the bosonic and fermionic 
matrices as a vector and as a Majorana-Weyl spinor, respectively. The bosonic matrices 
represent the dynamically generated space-time. A <i-dimensional space-time is described by 
configurations with d bosonic matrices having much broader eigenvalue distributions than 
the other (10 — d) matrices, up to some SO(10) transformation.^ Our four- dimensional space- 
time may be accounted for, if the d — 4 configurations (in the above sense) dominate the 
integration over the bosonic matrices. It was found recently that the IIB matrix model is 
indeed endowed with a natural mechanism that may realize such a scenario [^. [|. 

The realization of our space-time as a 'brane' in a higher-dimensional space-time has 
attracted much attention as an alternative to the more conventional approach in string 
theory using compactification (see Ref. [J7| and references therein). It turned out that such 
a set-up has many phenomenological advantages, including possible mechanisms which may 
solve the cosmological problem and the hierarchy problem. However, the dynamical origin 
of the brane has not been discussed so far. The IIB matrix model enables us to investigate 
whether a four- dimensional space-time emerges dynamically as a brane in ten-dimensional 

1 Dimensionally reduced Yang-Mills theories were first studied in Refs. Q. 

2 When one defines the IIB matrix model nonperturbatively, a Wick rotation to Euclidean signature is 
needed. Hence by Lorentz symmetry we actually mean SO(10) symmetry. 

3 Throughout this paper, we denote the initial dimensionality by D and the dimension after a possible 
SSB as d. 
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type IIB superstring theory through some nonperturbative effects Q. 

The spontaneous breakdown of Lorentz symmetry in matrix models has been addressed 
first in the bosonic case ||, where fermionic matrices are omitted (for recent work on the 



bosonic model, see Ref. |10|, |IT]). There the absence of SSB has been established by both an 
analytical method (to all orders in a 1/D expansion) and by Monte Carlo simulations. The 
same numerical result was obtained in the 6D and 10D SUSY matrix models |T2j , albeit with 
some simplifications to enable simulations at large N. The fermion integrals are complex in 
general in these cases, and the simulations were carried out including only the modulus, but 
omitting the phase. In addition, a low-energy effective theory was used in order to further 
reduce the computational effort. 

In Ref. |L3| we presented Monte Carlo simulations of the 4D version of the IIB matrix 
model, which is a supersymmetric matrix model obtained from the zero-volume limit of 4D 
super SU(iV) Yang-Mills theory. We were able to study the model with N = 16,24,32,48 
without any simplifications. These values of N turned out to be sufficiently large to extract 
the large N behavior of the space-time structure and to reveal the large N scaling for a 
number of Wilson loop correlators. 

Recently, it has been reported for the 4D SUSY model up to iV = 8 that the space-time 
is observed to be one-dimensional, if one selects configurations with large extent from the 
ensemble [|TJJ]. In the D-dimensional SUSY models in general, D = 4,6,10, configurations 



with large extent are suppressed only by the power — (2D — 5), independent of N [ 15 



Therefore, the observed anisotropic configurations may play some role in the large N limit, 
and such effects may also be relevant in other SUSY models, including the IIB matrix model. 

In this paper, we reconsider the issue of SSB of Lorentz invariance in the 4D SUSY 
model. If we adopt the conventional criterion based on the moment of inertia tensor, then 
the space-time appears one-dimensional, as suggested by the observation in Ref. flTJ|]. This 
would mean that the SSB does occur at large N. However, this conclusion depends on the 
definition of the order parameter, as we shall see. Thus we have to address the question 
which criterion for the SSB of Lorentz symmetry is actually physical. 

We recall that in the interpretation of the IIB matrix model as a string theory, the Wilson 
loops are identified with the string creation operators (T^|. Physical observables — such as 
scattering amplitudes — should be extracted from correlation functions of Wilson loops, 



4 A possible obstacle may be that gravitons propagate in ten-dimensional space-time, and hence one fails 
to reproduce the observed four-dimensional Newton's law. In Ref. || , it was demonstrated that this obstacle 
can be avoided in the case of D3-brane backgrounds due to the mechanism of Ref. M . 
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which were observed to have well-defined large N limits in D = 4 We therefore propose 
a "physical criterion" of SSB using Wilson loops as a probe, and we study it by Monte Carlo 
simulations of the full 4D SU(N) SUSY model. As a result, we find no trend of SSB up to 
N = 48. 



2 The 4D IIB matrix model 

The model we investigate is a supersymmetric matrix model obtained from the zero- volume 
limit of 4D SU(iV) super Yang-Mills theory. Its partition function is given by 

Z = J dA e~ Sb J dipdip e~ 5 ' , 
1 
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S b = -— 2 Tr[A„A 

S f = -^Tr^r^A^]) , (2.1) 

where A^ (fi = 1,...,4) are bosonic traceless N x N Hermitian matrices, and ip a , ip a 
(a = 1, 2) are fermionic traceless N x N complex matrices. The 2x2 unitary matrices 
are gamma matrices after Weyl projection; they can be given for example by 

T 1 = ia 1 , T 2 = ia 2 , T 3 = ia 3 , T 4 = 1 . (2.2) 

This model is invariant under 4D Lorentz transformations, where A^ transforms as a vector 
and ij} a as a Weyl spinor. The model is manifestly supersymmetric, and it also has a SU(iV) 
symmetry 

A^VA^ ■ ^ Q ^W f , $ a ^V$ a Vi, (2.3) 

where V G SU(iV). All these symmetries are inherited from the super Yang-Mills theory 
before the zero- volume limit. The model can be regarded as the four-dimensional counterpart 
of the IIB matrix model. 

The model is well-defined for arbitrary N > 2 without any cutoff. This was first con- 



jectured based on numerical results at small N confirmed further at larger iV [13| and 



and finally proved by Ref. |TI|. Therefore, the parameter g — which is the only parameter 



of the model — can be absorbed by rescaling the variables, 

A^g^X, ■ ij a =g 3 ^ a . (2.4) 

Therefore, g is a scale parameter rather than a coupling constant, i.e. the g dependence 
of physical quantities is completely determined on dimensional grounds. The parameter g 
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should be tuned appropriately as one sends N to infinity, so that each correlation function of 



Wilson loops has a finite large N limit. This issue has been studied numerically in Ref. [ |13|| . 
The conclusion is that the product g 2 N has to be kept constant when taking the large N limit. 



The tuning of g was also discussed in terms of analytical arguments in Refs. |Tj], [18], [19], ^(J . 



The integration over fermionic variables can be done explicitly and the result is given 
by det Ai, M. being a 2{N 2 — 1) x 2{N 2 — 1) complex matrix which depends on A^. The 
matrix M. is given by 

M aaW = (r M ) a/3 Tr (t a [A^t b }) , (2.5) 

where t a are generators of SU(iV), and we consider (aa) resp. (b(3) as one index. Hence the 
system we want to simulate can be written in terms of bosonic variables, 

-s b 



dAe~ bb detM. (2.6) 

A crucial point is that the determinant det M. is real positive [0 . This property was 
demonstrated in Ref. [ 131 , and it had been suspected earlier fPTf . (It also holds in other 4D 



SUSY models, see Ref. [[yj an d second Ref. in [Q.) Due to this property, we can simulate the 
model using a standard algorithm for dealing with dynamical fermions (the so-called Hybrid 
R algorithm [22]). In the framework of this algorithm, each update of a configuration is 



done by solving a Hamiltonian equation for a fixed "time" r. This algorithm is plagued by a 
systematic error due to the discretization of r that we used to solve the equation numerically. 
We performed simulations at three different values of the "time step" At and we extrapolate 
to At = 0. 



3 SSB of Lorentz symmetry ? 

In the IIB matrix model, the eigenvalues of the bosonic matrices A^ are interpreted as the 
space-time coordinates |2|, [18], [2~3| . We adopt this point of view in the 4D model as well. 
Since the matrices A^ are not simultaneously diagonalizable in general, the space-time is not 
classical. To extract the space-time structure we first define the space-time uncertainty A 
by 

^^^-^.fEIW)..} 2 . (3-D 

i 

which is invariant under Lorentz transformations and under the SU(A^) transformations 
( p.3| ). Formula ( |3.1| ) has been derived in Ref. based on the analogy to quantum mechanics, 
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considering as an operator acting on a space of states. As a natural property, A 2 vanishes 
if and only if the matrices are diagonalizable simultaneously. For each configuration A^ 
generated by a Monte Carlo simulation, we maximize A-iJJ )a\ 2 with respect to the 

SU(iV) matrix U. We denote the matrix which yields the maximum as ?7 max , and we define 

x in = (U max A^U^ nax )ii (3-2) 

as the space-time coordinates of N points Xi (i — 1, • • - , N) in four- dimensional space-time. 
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Figure 1: The eigenvalues of the tensor of inertia T as functions of the algorithmic "time 
step" At at N = 16, 24 and 32. The dashed lines on top represent the fits to eq. (\3.4) , 
confirming that the largest eigenvalue diverges logarithmically in the limit At — > 0. (The 
other lines are drawn to guide the eye.) 



In order to search the spontaneous breakdown of Lorentz symmetry, we first consider the 
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moment of inertia tensor of iV points x\. It can be defined as 

2 

Tfj,i/ — jY^jy _ ^ ~ %jfi)( x iu ~ x jv) 5 (3-3) 

which is a D x D real symmetric matrix.f] The D eigenvalues Ai > A2 > . . . > Xd > of 
the tensor T represent the principal moments of inertia. We take the average (A M ) over all 
configurations generated by the Monte Carlo simulation. 

In Figure [1] we plot the results for (X^/g against At at iV = 16, 24 and 32. (In all the 
Figures in this Section, we plot dimensionless quantities, so that they do not depend on the 
choice of the scale parameter g.) As At vanishes, (Xi)/g is observed to diverge as 

(Ai)/</~ - Cl lnAr + co (3.4) 

(where cq, c\ are constant in At), while the other eigenvalues converge.0 We recall that the 
extent of the space-time R defined by 

was found to be divergent [P~3|| . Since R 2 = ^AX^), at least (Ai) has to diverge as At — > 0. 
What we observe here is that (Ai) is indeed the only divergent eigenvalue. This is consistent 
with the observation in Ref. [14j . 



Let us introduce the probability distribution for the distance of two space-time points as 



f*r) = N{N iA >Jl>--\/^-->:i) 2 ) ) ■ «3.o 

Then R 2 can be written as 

poo 

R 2 = dr r 2 p(r) . (3.7) 
Jo 

The observed logarithmic divergence of R 2 is consistent with the asymptotic behavior 

p(r) ~ r" 3 , (3.8) 



5 This quantity has also been studied in Refs. 1191, p^j . Alternatively, one may define the moment of inertia 
tensor by 1^ = -hTr (A^A^), as it was done in Refs. || We have also measured the eigenvalues of that 
tensor, and the spectra are in qualitative agreement with the eigenvalues that we present here (based on 



definition (3.3)). 

6 We note that based on our data, a logarithmic divergence of (A2) is unlikely, but it cannot be absolutely 
excluded. This slight uncertainty is still a little stronger if we use the definition = -^Tt^A^A^) for the 
moment of inertia tensor. 
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Figure 2: The eigenvalues of the new tensor T^ w ^ as functions of the algorithmic "time 
step" At at N = 16, 24 and 32. The straight lines are the fits to eq. ($.1J\) . We see that 
now all eigenvalues converge in the limit At — > 0. 



which was predicted analytically fI5| . Based on this observation, a modified definition for 



the extent of the space-time has been introduced in Ref. |13 

2 



N(N — 1) 



i<3 



dr rp(r) 



(3.9) 



which turned out to be finite — as expected from relation fl3.8| ). The large N behavior of 
this quantity has been observed to amount to -Rnew/y^ = 3.30(1) • iV 1//4 , which is consistent 
with the prediction based on the low-energy effective theory 
This motivates us to define analogously a new tensor 



fii(iiew) 



(3.10) 
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Figure 3: The eigenvalues of the new tensor T^™ , extrapolated to At = and normalized 
by the factor (y/giV 1 / 4 )" 1 . We also include the results for N = 48, which are obtained at 
At = 0.002 (our statistics at other values of At is insufficient for a sensible extrapolation 
in this case). 

Let us denote the D eigenvalues of the tensor y( new ) as ;\( ncw ) > \ 

>■■■> A ( n ncw) > 0. 



(new) 



(new) 



2 ^ ^ /y D 

, all the eigenvalues are expected to converge. In 



Due to the relation ^ A 

Figure |2| we show the results for (^ cw> ) / y^g again at iV = 16, 24 and 32. We carried out 



an extrapolation to Ar = assuming the observable (9 (At) at small At to behave as (13| 



Q{At) - 0(At = 0) oc (At) 2 • | In Ar| 



(3.11) 



In Figure § we plot the extrapolated and normalized eigenvalues (A/T™ ) / (^fgN 1 ^) against 
1/N. (This is the normalization needed for a finite large N limit.) We observe that they 
move closer together as iV increases. Therefore we cannot recognize any trend for SSB. 



4 A physical criterion in terms of Wilson loops 

The results in the previous Section reveal a subtlety in the issue of SSB in the 4D SUSY 
model. The crucial question is whether any signal of SSB can be probed by physical quan- 
tities, such as scattering amplitudes. Therefore we have to reconsider how the IIB matrix 



model is interpreted as a string theory. In Ref. [0 it has been demonstrated that Wilson 
loops in matrix models can be identified with string creation operators in string theory. 
Hence Wilson loop correlation functions are the only objects with a direct physical inter- 
pretation in string theory. So we should ask whether any signal of SSB can be probed by 
Wilson loops. 

We recall that the extent of the space-time can be probed by the vacuum expectation 
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value (VEV) of the "Polyakov line" |EJ 



Pip) = ^Tr exp(ijy4„) , (4.1) 



where p M represents the total dimensionful (and hence "physical") momentum carried by the 
string.[| The VEV (P(p)) depends only on p = f y/p~p~ due to the SO (4) invariance. It starts 
at 1 for p = and drops down to zero at some value p. Then 1/p is a measure for the extent 



of the space-time. In Ref. | 13| , it was shown that the one-point function (P(p)}, as well as 
other Wilson loop correlators, converge to a certain function of p in the large N limit, if 
the scale parameter g is taken to be proportional to 1/ \/N. This means in particular that 
p oc {y/gN l / 4 )~ l . (If we set g = (V/48) -1 / 2 , as we do in Figures £| and [|, we find p ~ 0.7 
for N = 16, . . . ,48.) Therefore, the large N behavior of the space-time extent probed by 
the Polyakov line is 1/p oc ^/giV 1 / 4 , which is consistent with the result obtained from -R new 
defined in eq. ( [3.9|) . 

Let us formulate the SSB of Lorentz symmetry by using the Wilson loops as a probe. 
For each configuration we perform a SO (4) transformation = A^ U A U so that = 
jjTr(A^A u ) becomes diagonal : I = diag(o»i, u>2, <^3, W4), where uj\ > 0J2 > uj 3 > cj 4 > 0. 
Then we define 

P>) = ^Trexp(zp4) . (4.2) 

As a consequence of the diagonalization, the VEV's (P^(p)) (p, — 1, ... ,4) do depend on fi 
for finite N. If these functions are different even in the large iV limit, we may conclude that 
the SSB of SO (4) symmetry occurs. 

This method is analogous to the Ising model, where the magnetization (|M|) (M is the 
sum of spins divided by their number), serves as an order parameter for the SSB of Z2 
symmetry. Taking the absolute value of M corresponds to making an appropriate SO (4) 
transformation from A^ to A^. Note, on the other hand, that (M) = for any finite lattice 
volume even if the SSB takes place (for infinite lattice volume). Similarly, the VEV (P(p)) in 
eq. (fOD depends only on p = y/p^p^, even if the SO (4) symmetry is spontaneously broken. 

Our results are shown in Figure ^. (In this case we just present the results obtained at 

At = 0.002, which appears to be sufficiently small.) Here we set g = (iV/48)~ 1//2 and plot 

(P (U (p)) against iV at three different values of the momentum p, all of them below p. We 

7 This momentum p is identical to fc p h ys (= k/y/g) introduced in Ref. |l3|]. We take this opportunity to 
correct a typo in Fig. 5 of Ref. 1 13 . The label for the horizontal axis should be k/^/g instead of k 2 /g. This 
typo has propagated to our review articles j24| H|| as well. Fig. 3 of |Q and Fig. 3 of [B5| should have p 
instead of p 2 as the label. 



9 



observe no trend of SSB up to N = 48. To illustrate this observation in yet another way, we 
show in Figure [| the functions (P^(p)) for N = 16 and for iV = 48. We see that the results 
for the four Polyakov lines move closer together as N increases. 
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Figure 4: The four Polyakov lines (P^p)) (fi = 1, . . . ,4) atp = 0.316, 0.447, 0.548, where we 
take the scaling parameter g to be g = (A^/48) -1 / 2 . The Polyakov lines become approximately 
equidistant and they move closer together as N increases, hence we do not see any signal for 
SSB of Lorentz symmetry. (The lines are drawn to guide the eye.) 



Finally, we clarify the relation of the above result to the previous result ( |3.4j ). Let us 



denote the eigenvalues of A„ as a 



in 



of the a^i as 



Ux) 



\ i=i 



N) and introduce the probability distribution 
x - a^i) ) . (4.3) 
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Figure 5: The Polyakov lines (JF^(p)) (fi — 1, ... ,4) at N — 16 and at N = 48. T/je?/ move 
closer together as N increases, in agreement with the proposed scenario that they coincide in 
the large N limit. (The curves were measured quasi- continuously, and we show a few error 
bars as examples.) 



Then we can write (P^p)) as 




dxf^(x)e 



i px 



(4.4) 



Thus the VEV of the Polyakov line is just the Fourier transform of the distribution f^x). 
Our observation that (P^(p)) for fi = 1,2,3,4 approach the same function of p at N = oo 
implies that their inverse Fourier transforms / M (x) for /i = 1,2,3,4 also approach to the 
same function of x at iV — > oo. On the other hand, our result in the Figure [3] suggests that 
fi(x) ~ C(iV)/|a;| 3 at large x, but the / M (a;) with fi = 2,3,4 have a sharper fall-off at large 
x. This suggests that the coefficient C(N) vanishes in the large iV limit (at fixed g 2 N). 

We can rephrase the above statement in the following way. Note that (oj^) = (jjTr^A^) 2 ) 
(here we do not sum over the index fi on the r.h.s.) can be written as 




&xx 2 f^(x) 



(4.5) 



p=0 



The logarithmic divergence of (a>i) found in the previous Section (consider footnote [5]) implies 
that (Pi(p)) has a non-analyticity at p = of the form 



(Px(p)) ~l + C(N)p 2 \np + 



(4.6) 
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On the other hand, (P M (p)) with [i = 2, 3, 4 do not have this non-analyticity, since (002), (^3), 
(co> 4 ) are finite. Our observation that the four Polyakov lines (P^(p)) (/1 = 1, . . . , 4) converge 
to a single function of p suggests that the coefficient C (N) in eq. (|4.6|) vanishes in the large 
N limit (at g 2 N fixed). 

These statements can be checked by extracting the coefficient of the logarithmic di- 
vergence in eq. ( ^.4[ ). If we fit the data for (Ai)/g to the expected asymptotic function 
—Ci(N) In At + co(N), we find that the normalized coefficient C\(N) / \/ r N amounts to 2.2(3), 
2.52(7), 1.75(7) for N = 16, 24 and 32, respectively. Hence the behavior C,C — > as 
iV — > 00 is conceivable. 

To summarize, our observation in this Section implies that the quantities 

(/i = l,-..,4) (4.7) 

p=0 



f°° r 1 d r 

/ dxx 2 lim = - — lim (P^(p)) 



are all finite and equal. Notice that in eq. ( f4.7| ), the limit N — > 00 is taken before setting 
p = 0. The point is that one should first take the large N limit of the Wilson loop to make 
it actually physical, and then its derivatives at p = inherit a physical meaning, too. This 
does not need to be true for derivatives at p = for finite N. In fact, our results suggest 
that the p — > limit and the large N limit do not commute. 

5 Discussion 

In this paper, we wanted to clarify the issue of spontaneous Lorentz symmetry breakdown 
in supersymmetric matrix models, which was raised by Refs. |TJ]]. We propose a physical 
criterion for SSB, which we consider as a solution to this problem. In the particular case 
of the 4D SUSY model, configurations with only one- dimensional extent dominate when 
one adopts a conventional criterion for the SSB using the moment of inertia tensor, as was 
suggested by Refs. il4| . However, contributions of those configurations to physical quantities 



such as Wilson loop correlators seem to be strongly suppressed in the large iV limit. Indeed, 
if we rely on our physical criterion using Wilson loops as a probe, we do not observe any 
trend of SSB; the space-time probed by Wilson loops appears to be four- dimensional. 

Let us comment on the cases D = 6 and D = 10. We recall that in the D-dimensional 
SUSY models in general, D = 4, 6, 10, the eigenvalue distribution of has a slow fall-off 



with the power —(2D — 5) independent of iV [15|. Therefore, as pointed out also in Ref. [|l 



a problem may arise if one considers a tensor such as (I 71 )^, where = ^^(A^Ay) and 
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n > D — 3. The VEVs of the corresponding eigenvalues can be written as 

/ 1 N \ r°° H 2n 
<M n > = ( nJ2M 2h ) = J Jx**Mx) = (-if ^stA(p)> 



(5i 

p=0 



However, one has to take the large N limit before setting p = 0, in order to have a correct 
interpretation in string theory. The possible divergence of ((ui) n ) for finite N does not 
immediately imply an appearance of the SSB of Lorentz symmetry in physical quantities, as 
we have seen in D = 4. On the other hand, the conventional moment of inertia tensor 1^ 
(or T^ v in Section |3|) does not have such a problem in D = 6, 10, since the p — > limit and 
the large N limit commute. If one observes an SSB with the conventional moment of inertia 
tensor, it immediately implies an SSB in physical quantities. 

The absence of SSB in 4D SUSY model (with the physical criterion) is consistent with 
the conjecture that the phase of the determinant in the fermionic partition function plays 
a crucial role in a possible SSB of Lorentz symmetry. This conjecture is supported by the 
following results: 

1. The bosonic model does not show SSB ||. 

2. The SUSY model in D = 6 and D = 10 (using a low-energy effective theory) does not 
exhibit SSB if one omits the phase of the determinant |T3 . 



3. SSB does appear in the i/-deformed SUSY model in D = 6 and D = 10 in the large v 
limit, where v couples to the phase of the fermionic partition function ||. 

4. No SSB occurs in the 4D SUSY model (where the fermionic partition function is real 
positive), as discussed in the present work, where we do not use any simplification of 
the full model. 

When one integrates over the bosonic matrices in the 6D or 10D IIB matrix model, the 
phase of the fermion integral fluctuates rapidly except for the vicinity of configurations, for 
which the phase becomes stationary. In fact, this happens for any <i-dimensional configu- 
rations with 3 < d < D — 2 0. Those configurations are therefore considerably enhanced 
compared to the case where the phase is omitted. By using a saddle-point approximation, 
it was found that only the dimensionalities within the above range are possible. Note in 
particular that d — 4 is not excluded. 
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